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Spherical spin-orientation degeneracy of basic antiferromagnetic configurations 
due to the dipolar interaction in cubic lattices 

Eugene V. KholopcnQ 

Institute of Inorganic Chemistry of the Siberian Branch of the Russian Academy of Sciences, 630090 Novosibirsk, Russia 

Based on a simple general relation for the Lorentz field, the precise values are obtained for the 
energies of ferromagnetic and basic antiferromagnetic states in sc, bcc, fee, and diamond cubic 
lattices. Within both the 'nearest-neighbor' and group-theoretic approaches, a large variety of anti- 
ferromagnetic states revealed as corresponding to the lowest energies is addressed, with recognizing 
their spherical spin-orientation degeneracy, which turns out to be rather special in fee lattices. 

PACS numbers: 75.10.Hk, 75.25.+Z, 81.40.Rs 



The dipolar interaction of spins is typical of dilute mag- 
nets and so remains actual upon describing modern 
complex compounds where the direct exchange is either 
depressed 0,M 0, 0| of frustrated @. In this connec- 
tion, dipolar solutions in cubic structures are interest- 
ing themselves 0, S S and also enable one to un- 
derstand magnetic ordering in structures of a lower sym- 
metry [lCJ, [llj . It is known that antiferromagnetic 
ground states are typical of the sc and diamond |l3| 
structures, whereas the bcc and fee lattices are ferromag- 
netic HlHEi in thc thermodynamic limit [HlUlHiq. 

In order to describe basic spin configurations, the 
'nearest-neighbor' 01 an d group-theoretic Q treat- 
ments are usually applied 5]. However, an explicit motif 
of spin arrangement in shortest-distance directions [l2j |. 
that is characteristic of the dipolar interaction, is not 
still recognized in the bcc lattice. Nevertheless, the same 
motif gives rise to the widespread standpoint that basic 
antiferromagnetic arr ang ements are attached to crystal- 
lographic axes [3, LUJ Ljj • This is at odds with spherical 
degeneracy of ferromagnetic configurations Q , as well as 
with the fact that the spherical degeneracy of antiferro- 
magnetic states in sc lattices has already been proved 
[l4|. This is the reason that here we study this question 
addressed to all cubic lattices in a systematic manner, 
with treating spins as classical vectors. 

As for the problem of summation of dipolar lattice se- 
ries that is rather tedious conventionally 0, 0, 0] , here 
we propose a simple general approach based on intro- 
ducing fictitious charges . Let us consider the general 
case of a triclinic crystal with lattice constants a, b, and 
c along non-orthogonal triclinic directions e a , e b , and e c , 
respectively, and with n moments Mj localized at the 
positions tj in a unit cell, providing that Mj can differ 
from one other. Concentrating on a jth sublattice built 
up of parallel moments Mj at a given j, it is expedient 
to cast Mj in terms of triclinic coordinates 

Mj = M io e + M jb e b + M jc e c = [M ja , M jb , M jc ] , (1) 

where the square brackets are used in order to distinguish 
this representation from a Cartesian one. Then the unit 
cell of this sublattice, with Mj at its center, can be mod- 



ified by three point charge species 

Mia 



qji 



b ' 



Qj3 



M ic , . 
(2) 



at the unit-cell positions 



r 1 = [ia,0,0], r 2 = [0,|6,0], r 3 = [0, 0, §c], (3) 

but the charges — <2ji, ™<7j2j an d — g j3 are assumed at 
— ri, — r2, and — r3, respectively. As a result, one can 
see that the total dipolar moment of the unit cell is zero 
and all the charges above cancel one other upon sum- 
ming over the lattice. Therefore, the field exerted by this 
modified sublattice on any moment M at the position 
t is determined by the convergent lattice series, which 
can be easily obtained from the appropriate expansion of 
a Coulomb lattice series, as will be discussed elsewhere 
[2l| . This is a generalized Lorentz field of the form 



i 



3((R,-t)M J )(R,-t) 



M, 



|R,-t| 5 



|R< 
t 



1*777, ^ 



IRi 



},(4) 



where the summation is over the jth sublattice described 
by Ri, the prime on the summation sign implies missing 
the cases of zero denominators, ((R, — t)Mj) stands for 
the scalar product. The first term in the curly braces 
in (0} is the direct sum of dipole contributions and the 
second one is the Lorentz field itself that is responsible 
for the topological boundary effects Note that the 
convergence of (@J takes place if the summation over i 
is the outer procedure. Keeping Q in mind and taking 
all the sublattices into account, the total bulk dipolar- 
energy per unit cell is as follows: 



£ = -i^(M il Lj 3 (tj 1 -tj 2 )). 



(•5) 



Here we focus on the four cubic structures with the lat- 
tice spacing a: the sc, bcc, fee, and diamond ones. The 
volume v s per spin is equal to a 3 , a 3 /2, a 3 /4, and a 3 /8, 
respectively. All the spin values |M| are supposed to be 



2 



TABLE I: The bulk dipolar energy £ s per spin in units of 
M. 2 /v s for the ferro- (F) and basic antiferromagnetic (AF) 
states in cubic structures, with indicating the three- (3D) and 
two-dimensional (2D) degeneracy in the latter case. 



Structure 


£, (F)° 


£s (AF) 


AF dc 
3D" 


generacy 
2D" 


sc 


-2.094395 c 


-2.676789 


l d 




bcc 


-2.094395 c 


-1.985921 c 


2 


3 d 


fee 


-2.094395" 


-1.807574 


4' 


3 d9 + 3' 1 


diamond 


-2.094395 c 


-3.196851 


2 


3 d 



"The direction of a reference spin species is arbitrary. 

'Reference spins lie in one of the three basal planes. 

This is just the value of —2n/3, in accord with JHJ 14]. 
H ere t he c ases are contained of co lline ar spins inherent in the sc 
|l2l . bcc I'ill . fee llij . and diamond lattices. 

e This value was discussed earlier Il4l as associated with spins 
along different crystallographic axes in the two sc sublattices. 

^Two-parameter sets of states at a given reference spin. 

9 One-parameter continuous sets of states. 

h One-parameter piecewise continuous sets of states. 



the same, so that the energy per spin £ s = £ jn is of in- 
terest. The energies for ferromagnetic spin arrangements 
calculated from J2J)-© are shown in Table [I] and agree 
with the fact that the direct dipole contribution to Q 
upon summation over cubic shells is zero by symmetry, 
so that the bulk dipolar energy per spin is reduced to the 
contribution of the Lorentz field itself |22j: 



£ s = 



2ttM 2 

3w, 



(6) 



This approach is also fruitful upon calculating basic 
antiferromagnetic states, which can be decomposed into 
sc ferromagnetic sublattices. The results are exhibited 
in Table U as well. To address a large variety of such 
states revealed, we resort to the description of Luttinger 
and Tisza |l4| , in terms of which the sc antiferromagnetic 
structure is described by the state 

(aX 5 + 0Y 5 + 7 Z 5 ), a 2 + /3 2 + 7 2 = 1 (7) 

with the energy per spin proportional to |l4| 

a 2 +/3 2 + 7 2 . (8) 

According to the definition of X5, I5, and Z5, it implies 
the three-dimensional (3D) degeneracy of this antiferro- 
magnetic (AF) state shown in Fig. |T1 (a), with local mo- 
ments, in units of |M|, of the form [F4J 



a = (a,/3,7), 
c = (-a,/3, -7), 



b= (0,-/3,-7), 
d = (-a,-/3, 7). 



(9) 



The basic AF bcc structure is unstable, in agreement 
with Table |U Nevertheless, it can be specified in the 
same manner. If the two sc sublattices are described by 
the states 



> 



a c 
b 



a (3D) 



"c e" 

bl (3D) 



,a n 
-n< 



K n 



b2 (2D) 





cl (3D) 



7\ — , A 
■af fa 




c2 (2D) 
-a 



d (3D.2D) 



FIG. 1: 3D basic antiferromagnetic configurations in the sc 
(a), bcc (bl), and diamond (cl) structures, with spins speci- 
fied by © and (1121) . 2D configurations typical of the bcc (b2) 
and diamond (c2) structures contain spins driven by 11411 . A 
basic antiferromagnetic fee structure (d) is common to 3D and 
2D configurations described by I15H and either by 1201 - 12311 
or by I|28|l - lj30|l as an example, respectively. 



normalized like Q , then the energy of interaction of these 
sublattices is specified by 

max{«77 + /3C + 7C} ~> £ = 7, V = a, C = P, (H) 

for the expression in the curly braces can be treated as a 
scalar product of two unit vectors which must be paral- 
lel. It implies the 3D configuration shown in Fig.nj(bl), 
where apart from @, the local spin states defined as 



e=(7,a,/3), f = 

g= (-7, a, -j3), h 



(j,-a,-(3), 
: (-7,-a,/3) (12) 



(aX 7 + (3Y 7 + jZ 7 ), (£.X 6 + V Y 6 + (Z 6 ) 



(10) 



take place. The interchange of the two sublattices leads 
to the twofold degeneracy of such states, as shown in Ta- 
ble 0] In other words, the two rotations of the lattice in 
Fig. n ( a ) about the (1,1,1) axis, but without rotating 
spins, generate the first sublattice here. It is worth not- 
ing that the particular cases of a along one of the four 
cubic diagonals are associated with the idea of Sauer 01 



3 




FIG. 2: The structure of interpenetrating ferromagnetic 
chains shown by the solid arrows in the nearest-neighbor di- 
rections as a particular case of Fig.0(bl) at a = (-±= , -±= , -±=) . 
The first sublattice is depicted by the dotted lines, whereas 
the filled circles exhibit the second sublattice. 



about ferromagnetic chains along the nearest neighbors, 
as shown in Fig. [51 where the cubic chain packing ap- 
pears to be similar to the chain cubic structure intrinsic 
in TIHF2 at low temperatures [2^ |. 

Another possibility is associated with the states 

(/3Y 7 + lZ 6 ), ( 7 Y 7 + f3Z 6 ), (13) 

which imply the configuration shown in Fig. ^ (b2) at 



a = (0,/3, 7 ), n=(0,7,/3). 



(14) 



The cyclic interchange of the parameters in (|14|l explains 
the threefold degeneracy of this two-dimensional (2D) 
set, as indicated in TableQ] Note that all the states above 
follow from the 'nearest-neighbor' treatment of Sauer 0] 
as well and this set of possibilities is complete. 

The same 'nearest-neighbor' treatment enables us to 
extend the solution typical of the bcc structure to the 
diamond one. Indeed, any sc lattice can be regarded 
as a combination of two fee sublattices, as it happens 
in the NaCl structure, but with equal energy effects at 
body-centered positions. So, upon removing one of the 
fee sublattices from either of the sc sublattices in the bcc 
structure, a diamond structure arises. As a result, the 
total set of solutions for the diamond structure is still 
defined either by © and 1)120. or by (|14ll . as shown in 
Fig. n (cl) an d Fig. ^ (c2), respectively, with the cor- 
responding degrees of degeneracy. The chain structure 
similar to that in Fig. but with two alternating spin- 
spin distances along every chain, takes place here as well. 

For completeness, we also consider the fee lattice, 
where the situation turns out to be much more compli- 
cated and dramatic. In terms of unit vectors a, ni, 112, 
and ri3 specifying four sc sublattices of the same form, 
both the above treatments predict the limiting antiferro- 
magnetic structure drawn in Fig. (d) . On introducing 
the unit vector components along the principal axes: 



one can show that the lowest energy of such a structure 
corresponds to the condition 

(m + Ca - a) 2 + (& + C2 - P) 2 

+ (6 + ?72-7) 2 + (6+% + Ci) 2 = 0, (16) 



with the obvious solutions 

m + C3 - a = 0, 
6 + m - 7 = 0) 



6 + C2 - /3 = 0, 

6 + V3 + Ci = 0. 



(17) 



Further it is expedient to write down the unit vector nor- 
malization in the form 

6 = A, m = CiP h Ci = OiQi, (18) 



Af + Cf 



1. 



p 



(19) 



corresponding to spherical coordinates. According to 
()16J) and (|19|l . it is clear that apart from the reference 
unit vector a, which is assumed to be arbitrary, two ad- 
ditional parameters describing l|18fl remain variable as 
well. In particular, we choose £1 = A and £2 = B as such 
parameters here. On substituting 1(18(1 into (|1T|1 and re- 
solving those relations with account of 1(19(1 , one can show 
that the other components take the form 



C2 = P - A 



'/2 



= ±yi - (fi - a) 2 — b 2 



& = 7 T V~l - (/3 - A) 2 - B 2 ' 

-a{A 2 - S 2 - a 2 - B 2 ) ± Bfl 
2 (a 2 + B 2 ) ' 
B(A 2 ~£ 2 -a 2 -B 2 )±aQ 



(20) 
(21) 



Ci 
m 

Cs 



(22) 



2(a 2 + B 2 ) 
—B(A 2 - + a 2 + B 2 ) t aft 

2{a 2 + B 2 ) 
a(A 2 -gf + a 2 + B 2 )tBQ 
2{a 2 + B 2 ) 

Q = V /4(1-A 2 )(a 2 +-B 2 ) - [A 2 -£ 2 -a 2 -B 2 ) 2 , (23) 

where the choice of an upper or lower sign is common to 
both the relations (|21(l and the same is right for all four 
relations (|22|l . which case is independent of 1)21(1 . Alto- 
gether, four two-parameter sets of solutions as functions 
of a, /3, and 7 exist, as pointed out in Tabled The change 
in A and B is restricted by the following inequalities: 



\A\<1, 

{a 2 e 



T 2 7 - 7 2 < 09 - 
y 2 - B 2 ) 2 < 4(1 



A) 2 + B 2 
A 2 )(a 2 + 



< 1, 
B 2 ) 



(24) 
(25) 



(a, 0,7), n; = (CuVhCl), 



(15) 



which emerge from l(19|l and from the conditions that the 
square roots in l(21[) and l(23[) be real, an upper or lower 
sign in (|24|l agrees with that in 121|) . 

In particular, if A = f3 and B = A = ±1, then relations 
and (|5!)|l hold automatically and solutions (f^Uf) — (f?5f) 



4 



(2 = 0, & 



a(l - 


-/3 2 )±AV(1- 


/? 2 )(1 


-7 2 ) 




1 + a 2 






A(/? 2 


-i)±<V(i- 


/3 2 )(1 


-7 2 ) 




1 + a 2 






A(7 2 


-1)T"V(1- 


(3 2 )(l 


-7 2 ) 




1 + a 2 






a(l - 


7 2 )TAV(1- 


(3 2 )(1 


-7 2 ) 


1 + a 2 



(26) 



(27) 



are reduced to 

m = 
m = 
Ci = 

'?3 = 

c 3 = 



The existence of four branches of solutions l|20(l - (|23|l en- 
sues from by continuity. 

The case of a = B = implies a = (0, f3, 7) and K2 = 
(0, r)2, Ca)- It is degenerate in the representation at hand 
and can be studied directly from (|17l) giving rise to 



V3 



(3 



-Vi, 



Vl + Ci = 1 - & (28) 



According to i|19[l . there are two solutions for the other 
components here. The first one is defined by 



6=6 = 0, 



m = 7, (2 = P, 



(29) 



in which case the consistent rotation of a and H2 in the 
plane normal to the first axis is accompanied by an inde- 
pendent consistent rotation of the couple of ni and 113 in 
the same plane. In other words, 1|28|) and (|29|) generate a 
one-parameter continuous solution. The second solution 
in question is of the form 



6 = -6 = /? -7, n 2 = a at ^7 > 0, 
6 = 6 = P + 7, n 2 = -a at 07 < 0. 



(30) 



According to if^Sl) and (fSUf) , one can see that an indepen- 
dent consistent precession of ni and n 3 about the first 
axis is admissible, while both [3 and 7 are nonzero and 
so determine the angle of precession. In this case the 
rotations of both a and n.2 are connected. But if either 
[3 or 7 goes through zero, then H2 undergoes the sudden 
inversion. At those moments the precession of ni and n.3 
disappears: both ni and 113 are along the first axis and 
one of them is inverted in a jump-like manner as well. 

Note that the solutions of similar types are to be ex- 
pected at (3 — 773 = and 7 = d = 0. This is the 
reason that all three pairs of these one-parameter two- 
dimensional solutions arc pointed out in Table [IJ though 
the latter solutions, with including piecewise continuous 
ones, are contained in the general case of 120|) (|23ll . 

In summary, both the 'nearest-neighbor' and group- 
theoretic approaches are suitable for investigating the 
lowest energy antiferromagnetic configurations. The lat- 
ter one appears to be more effective upon describing the 



sc and bec lattices, but the former is favorable for extend- 
ing the foregoing results to the diamond case. Of course, 
the spherical degeneracy discussed is the most important 
for antiferromagnetic ground states in the sc and dia- 
mond lattices, where it accounts for the high isotropic 
susceptibility recorded experimentally 0] . On the other 
hand, the spherical degeneracy revealed is a general prop- 
erty of the cubic symmetry. Thus the peculiar features 
of the foregoing results relevant to the sc, bec, and fee 
lattices are to be important for understanding the energy 
splitting upon lowering the lattice symmetry. 
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